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We analyze theoretically the dynamics of degenerate condensate of cold indirect excitons. We 
account for both linear spin dependent terms arising from spin-orbit interaction of Rashba and 
Dresselhaus types and non-linear terms transforming a pair of bright excitons into a pair of dark 
ones. We show that both terms should lead to the qualitative changes in the dynamics of cold 
exciton droplets in the real space and time. 



I. INTRODUCTION 

Collective phenomena lie behind many remarkable ef- 
fects in physics. One of their famous manifestations is 
Bose- Einstein condensation-^ which occurs if a system of 
Bose particles is cooled down beyond the critical temper- 
ature Tc which strongly depends on properties of the in- 
dividual particles, in particular their effective mass. For 
systems of cold atoms, where condensing particles are 
very heavy, Tc lies in the nano-Kelvin regime, which rules 
out any possibility of the practical implementation of this 
phenomenon. 

On the other hand in the field of condensed mat- 
ter physics various candidates were proposed for the re- 
alisation of BEG with critical temperatures orders of 
magnitudes higher then those of cold atoms .'^ The for- 
mation of exciton condensates in bulk semiconductors 
was theoretically predicted more than 40 years agOj^^ 
but appeared to be difficult to realize experimentally. 
Since then, other solid-state candidates were proposed 
for achieving high-temperature BEG, including Quantum 
Hall bilayers,!^ magnonsP cavity exciton-polaritonS^'^ 
and indirect excitons .'^ESl xhe latter system is in focus 
of the present paper. 

Spatially indirect exciton is a bound state of an elec- 
tron and a hole localized in coupled parallel 2D layers. 
Electron and hole wave functions show a very small over- 
lap and consequently indirect excitons have a long life- 
time as compare to ordinary excitons. They behave as 
metastable particles which allows their cooling beyond 
the temperature of a quantum degeneracy ! ^ ^ I ^ ^ I 

The indirect excitons have been widely studied both 
experimentally and theoretically in recent years. Super- 
fluid behavior of a system of indirect excitons was pre- 
dicted by Lozovik and Yudson more than 30 years aj 
and subsequent theoreticaP^'ifi' and experiment aP^^^i^l 
studies have suggested that this should be manifested in a 
series of remarkable effects, including persistent currents, 
Josephson-related phenomena and spontaneous pattern 
formation in the real space. 

Surprisingly, most of the works dedicated to indirect 
excitons have neglected their spin structure. On the 



other hand, it became clear in recent years that the ac- 
count of the spinor nature of the condensing bosons can 
lead to new qualitative phenomena. For cavity polaritons 
accounting for the spin led to the appearan ce of spinop- 
tronics, an optical analogue of spintronics.'22l2ll It was 
also shown that spin dependence of polariton-polariton 
interactions can lead to the appearance of intriguing non- 
linear polarization phenomena in polariton condensates, 
such as polarization multistability,^^ full paramagnetic 
screening also known as spin Meissner effectPSl and spin- 
dependent condensate velocities in the hybrid Bose- Fermi 
systems. One can expect that the spinor structure of 
indirect excitons c an als o have dramatic impact on their 
collective behavior.f^SllSl 



II. THE MODEL 

The spin of an indirect exciton is inherited from spins 
of the individual electron and heavy hole forming it. 
The possible spin projections of the electron's spin on 
the structure growth axis (z-axis) are ±1/2, while possi- 
ble spin projections of heavy holes spin are ±3/2. The 
exciton thus can have four possible spin projections, 
= ±1,±2 (Fig. [l] (a)). The bright states with 
5^ = ±1 can be created by external right or left cir- 
cular polarized light, while optical creation of the states 
with Sz = ±2 is prohibited by selection rules. How- 
ever, these states known as dark states cannot be ex- 
cluded from the consideration, as they can appear due 
to the presence of spin-orbit interaction (SOI) of Rashba 
or Dresselhaus type or can be created as a result of col- 
lision of two bright excitons with opposite circular po- 
larizations, as will be discussed below. For direct ex- 
citons the energies of bright and dark states are split 
off by electron-hole exchange interaction with character- 
istic value of tens of microelectronvolts.— However, for 
indirect excitons, where the overlap between the wave 
functions of the electrons and holes is very small, this 
splitting can be neglected and bright and dark excitons 
can be considered to have the same energy. 

Below the temperature of the quantum degeneracy, a 
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system of cold indirect excitons can be thus described by 
a four component macroscopic wave function 4'(r,i) — 
('i'+2(r, t), ^'+i(r, t), *_i(r, i), *_2(r, t)), where the sub- 
scripts correspond to the z-projection of the spin. Its 
dynamics can be obtained from the following equation 



masses and mimics the renormalization of indirect exci- 
ton dispersions due to spin-orbit interaction of electron 
(see Appendix A for details). Going to the real space 
representation one can use the explicit expressions for 
momentum operators — —idx, ky — —idy. 



SH 



(1) 



where H represents the Hamiltonian density of the sys- 
tem accounting for free propagation of particles, SOI 
and exciton-exciton interactions and can be represented 
as sum of single-particle and interaction parts H ~ 
Hq + Hint ■ Lst us consider the terms Hq and Hint sepa- 
rately. 

The term Hq can be calculated as 



i/o = 4'^(r,t)t*(r,i), 



(2) 



where the 4x4 matrix T contains terms correspond- 
ing to the kinetic energy and interactions with effective 
magnetic fields of various types. The latter can be di- 
vided into three categories. First, there is SOI acting on 
the spin of the heavy hole. It scales as a cube of the 
kinetic momentunP^ and we neglect it in further con- 
sideration. As for the SOI acting on electrons spin, it 
can be represented as a sum of the Dresselhaus term 

Hd = fi {^a^kx — o'yky^ arising from the lack of inversion 

symmetry for the crystalline lattices of most common 
semiconductor materials (GaAs, CdTe etc) and Rashba 

term Hn = a (^<Txky — cTykx^ appearing due to the struc- 
tural asymmetry of the QW in z-direction.'22l The coefH- 
cients a and /3 are constants which depend on the mate- 
rial and geometry of the structure, and a^^y are Pauli ma- 
trices. For the indirect excitons both Rashba and Dres- 
selhaus terms lead to the transitions ±1 — ±2 mixing 
bright and dark exciton states and matrix T thus reads 



Ti2 

ti2 



(3) 



Now let us consider the term accounting for exciton- 
exciton i nterac tions. As indirect exciton is a compos- 
ite boson,'^2EII they can be divided into four categories, 
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FIG. 1: (Color online), (a) Spin structure of indirect excitons. 
(b) Feynmann diagrams showing different types of possible 
interactions between two excitons. 



where 2x2 blocks T12 can be written in momentum 
space representation as 



■ 12 



h^k^/2M 



(4) 



with the operator Sk 



X 



I3{kx -f iky) + a{ky + ikx) 



X 



me/M is the ratio of effective electron to exciton 



namely the terms corresponding to the direct Coulomb 
repulsion, exchange of electrons, exchange of holes and 
simultaneous exchange of electron and hole (exciton ex- 
change). These processes can be visualized using the in- 
teraction diagrams shown at Fig. [l] (b). 

The corresponding interaction Hamiltonian thus reads 



a=±l±2 

x{\^+^\''\^+2? + |*-lP|*-2|' + \^ + l\''\^-2? + \^-l\''\^+2?) + {Ve + K) (^l *- 1 *+2 2 + 4';2 *-2 * + l ) ,(5) 
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where Vdir, Vx, Ve and Vh denote direct dipole-dipole re- 
pulsion, whole exciton exchange, electron exchange and 
hole exchange Coulomb interaction, respectively. Since 
we are interested in the behavior of a weakly interacting 
Bose- Einstein condensate of indirect excitons, the main 
contribution comes from the processes with zero trans- 
ferred momentum, and all values of the matrix elements 
in the above expression are taken for g = 0. It is well- 
known that contrary to the case of conventional exci- 
tons, the spin-independent direct interaction of indirect 
excitons does not vanish for zero excha nged mo menta 
due to the strong dipole-dipole repulsion .MEIIIII in ad- 
dition, the processes of electron and hole exchange also 
have influence on dynamics of the system making it spin- 
dependent. 

The first term of the Hint corresponds to the first line 
of interactions in the Fig. [l] (b) and describes all pos- 
sible interactions between indirect excitons of the same 
spin configuration. The second term corresponds to the 
processes of the direct Coulomb interaction and exciton 
exchange shown at the second line. The third and fourth 
lines of interaction diagrams can be combined in the third 
term of Hint ■ Finally, the fifth line in the Fig. [l] (b) cor- 
responds to the fourth term of Hint and leads to the 
transition between pairs of bright and dark indirect ex- 
citons. 



The interaction constants corresponding to all four 
types of interaction can be estimated in the same fashion 
as for direct excitons. The direct Coulomb interaction 
between indirect exciton is given by 



Vdir — 



eeo ' 



(6) 



where L is the distance between centres of QWs. The es- 
timation of electron and hole exchange interaction con- 
stant is more cumbersome and requires numerical cal- 
culation of an exchange integral (B17) given in the Ap- 
pendix B (see also Ref. [33]). 

For indirect excitons in the long wavelength limit 
(g — > 0) direct and exciton exchange interactions coin- 
cide (Vdir = Vx) as well as electron and hole exchange 
(Vg — Vh)- It should be noted that the electron and 
hole exchange interaction strongly depends on the dis- 
tance between the centres of QWs and changes the sign 
for certain separation, as it is shown at Figj2]). This fact 
is important for further consideration. 

The set of spinor Gross-Pit aevskii equations describing 
the dynamics of the system and accounting for both lin- 
ear terms coming from Hq and non-linear terms provided 
by particle-particle interactions of various types can be 
obtained using Eqs. |l|5|and reads: 



dt 
dt 



ih 
ih 



= - + Vo^+i\^+i\^ + (Vb - w)-^+i\^-i\^ + yo*+i(l*-2|' + 1«'+2|') + w^*li*+2*-2, (7) 

= + 5i2«'-2 + Vo^-l\^-l\^ + {Vo - T4^)*_i|* + i|2 + Fo*_l(|*+2|' + |*-2|') + I^*;i*+2*-2, (8) 

= m+2 + ^12^+1 + Vb*+2|*+2|' + (Vb - T4^)^+2|*-2|' + yo*+2(|*-i|' + |*+l|') + VK^^a^+l*-!, (9) 

= ^*_2 - 5i*2*-l + K)^-2|*-2|' + {Vo - M^)^_2|*+2|' + Fo*_2(|*+l|' + + I^^';2*+l 1 (10) 



where Vb = Vdir + Vx + K + Vh, W = Ve + Vh and 
we defined the kinetic energy operator E — —h'^W^/2M. 
Due to the symmetry between electron and hole exchange 
at g = the set of the equations we consider contains 
only two independent interaction parameters Vq, W. This 
differs from the case considered in Ref. [25] where this 
symmetry was not accounted for and 5 parameters were 
used for description of the interactions in the system. 
In realistic system not all of these parameters will be 
independent. 



III. SPECTRUM OF ELEMENTARY 
EXCITATIONS 

The calculation of the spectrum of elementary excita- 
tions in the system can be done using the linearization of 



spinor Gross-Pitaevskii equations ([7|-([l0| with respect 
to small perturbation around the ground state. The spin 
configuration of the ground state of the condensate can 
be found by minimization of its free energy 

F(* + i,*_i,*+2,*_2,A*) = ^-M/(* + 1,*-1,*+2,*-2), 

(11) 

where fj, denotes the chemical potential of the condensate 
and/(4'+i,4'_i,^'+2,*_2) = |*+iP + |*-iP + |*+2p + 
|^'_2p. The condition /(*+i, *+2, 'I'-2) = n, 
where n is a total concentration, gives an additional equa- 
tion for the determination of /i. 

Let us investigate the generic Hamiltonian H — Hq + 
Hint more precisely. Its first part Hq consists of kinetic 
energy and spin-orbit coupling terms. Both of them de- 
pend on the velocity of the particles and can be usually 
disregarded when considering the ground state, where the 
interaction Hamiltonian (15]) plays a major role. Never- 
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FIG. 2: (Color online). Electron direct (blue dashed line) 
and exchange (green solid line) interaction for indirect ex- 
citons as a function of separation distance between quan- 
tum wells L. The sign of the exchange matrix element 
changes at a certain separation distance. We consider the 
case of GaAs/AlGaAs/GaAs heterostructure in the narrow 
QW limit. 



theless, we will later show that accounting for SOI leads 
to qualitative changes of the ground state. 

The symmetry of Hint with respect to the interactions 
between components results in a non-trivial ground state 
solution. Let us consider four particular cases of homo- 
geneous condensates: 

1) \^%\ - V^, *'Li,±2 = or \^^_,\ = ^, = 



or \^%\ = yn, = or \^%\ = ^, = 



(one-component condensate), 




2) = = or = 
'^j-i = ("m" two-component condensate), 

3) = v&o = or \^\^^\ = 
'^^i _2 — {"ij" two-component condensate). 



4) l^'j-i -1-2 1 ~ \AV4 (four-component condensate), 

where we defined the ground state wave function for 
each component. For instance, the first case implies the 
situation where only one spin component is present in 
the condensate with total density n, while the last case 
implies equal distribution of the density between the four 
spin components. 

Looking at the interaction Hamiltonian ^ one can 
note that the last spin-flip term is important. While all 
other terms do not depend explicitly on the phase of the 
condensate wave function, it does. Moreover, it is always 
possible to make the term 'if*^i'^*_i'i> +2^^ -2 negative by 
adjusting the relative phases of all four condensates, and 
thus it always lowers the energy of the ground state in- 
dependently of the sign of electron exchange term W. 
However, the sign of W will still affect the ground state 
properties since it enters into Vq matrix element. Thus, 



one should consider separately two different situations 
with negative and positive electron and hole exchange 
interactions. 

For negative W = Ve + Vh < (large QWs separations) 
the free energy of the four condensate states considered 
above is given by 

H(^) ^{V,„ + Vx + V, + Vh)'^ = ^, (12) 



H 



(2) 



{Vd^r +Vx + 



W 



2 ^)T-(^°-T)y' (13) 



= {Vd,r +Vx+Ve + Vh)'— = 



2 2 

2 ^r„2 



(14) 



H^'^ ^{Vd.r + Vx + V. + V^.)'^^^, (15) 



where we denoted by H^^\ Hlf, Hlf and Tf^^) the 
free energy of one-component, "m" two-component, "zj" 
two-component and four-component condensates, respec- 
tively. One can see that the ground state is 7-time 
degenerate and configurations of one-component, two- 
component and four-component condensates are possi- 
ble. The chemical potential in this case is equal to 
H< = [Vdir + Vx + Ve + Vh)n = Von. Therefore, the 
ground state of the system will be chosen by spontaneous 
symmetry breaking mechanism. The high level of the de- 
generacy of the ground states means that the system can 
demonstrate a large variety of the topological excitations 
(solitons, vortices and skyrmions). Moreover, in the sys- 
tem one can in principle observe the fragmentation of 
the condensate into domains with different spin struc- 
ture. The analysis of these interesting effects, however, 
lies beyond the scope of the present paper. 

For positive W = + Vh > exchange interaction 
(small QWs separations) the free energy of the system 
for different types of the condensates yields 



H(^)^(^Vd„ + Vx + V. + V,)'^^^, 



(16) 



-[.Vd^T + Vx^ ^ )^ - —) — , (17) 



Hlf = {Vd^r +Vx+Ve + Vh) 



2 ' 



(18) 



H^^^^iV.r + V. + ^f^-iVo~^4. (19) 

In this situation the ground state is three times degen- 
erate and ether four-component or two-component con- 
densate in -t-1, —1 or +2, —2 configurations is preferable. 

The spectrum of elementary excitations in the system 
can be calculated using the standard method of lineariz- 
ing the Gross-Pitaevskii equations ([7|-([l0l) with respect 
to small perturbation taken in form of a plane wave. For 
instance, in the four component condensate case with 
equal fraction of each spin state it is taken in the form 
= yrj74-hA,e^('""-'^*^+S*e-*('''-'^*).l2l]The solution 
of the system for small amplitudes Ai and Bi gives the 



5 



dispersion relations of the quasiparticles in the conden- 
sate. 

First, we consider the case when spin-orbit interaction 
is absent. For large separation between quantum wells 
the exchange interaction is attractive and the chemical 
potential of the ground state is defined as /i< = Von. 
The corresponding spectrum of excitations for interac- 
tion constants Vdir = 19.9 ^eV nm^ and 14 — —1.78 
fieV ^m? is plotted in the Fig. p[a). It contains linear 
Bogoliubov mode, gapped quadratic mode and double 
degenerated gapless quadratic mode given by relations 



Ek+n\Wl 



hjjJ2 



(20) 
(21) 
(22) 



where Ek = h'^P/2M. 

For the case of small separation between QWs the ex- 
change matrix element Ve changes sign and becomes re- 
pulsive (see Fig. [2]). The chemical potential of the ground 
state is thus defined as /j,> = (Vb — W/2)n. We plot the 
spectrum of elementary excitations for L = 6 nm sep- 
aration between wells where Vdir = 9.95 ^eV ^.rn? and 
Ve = 1.24 fj.eV ^m^ (Fig. ^c), solid lines). The disper- 
sion relations have the form 



^^i3A = \/Ek[Ek + 2/i>], 



huj> = ^Ek{Ek + nW), 



(23) 
(24) 



where the modes /i^i 2,4 are now three times degenerate. 
The obtained spectrum coincides qualitatively with the 
general dispersion relations given in the Ref. |25| . 

The account of the spin-orbit interaction requires in- 
cluding the Rashba and Dresselhaus terms in the Gross- 
Pitaevskii equations. One can expect that SOI removes 
the spin degeneracy and leads to the splitting of the 
modes. Moreover, due to the spin-orbital interaction one- 
component condensate ground state is no more prefer- 
able since the lowest energy state requires the presence 
of either (+l,+2) or (—1,-2) components. Here we 
present the numerical calculation of the quasiparticle 
spectrum accounting for the isotropic spin-orbit inter- 
action of Rashba type. For non-interacting case the dis- 
persion minimum moves to the fcg = ^^r- points which 
corresponds to the non-zero condensate phase velocity 
in the ground state (note that group velocity defined as 
Vg = hr^dE/dk remains equal to zero and condensate 
is not moving). The ground state of the condensate in 
this case acquires the total non- zero phase e**'"'', where 
the orientation of kg vector is defined by spontaneous 
symmetry breaking process. 

The corresponding renormalization of quasiparticle 
dispersion now occurs in the vicinity of kg points show- 
ing the linear spectrum (Fig. [sj^b, d)). This situation 
is reminiscent to the renormalization of the bogolon dis- 
persions in the exciton-polariton condensate where the 
role of spin-orbital interaction is played by longitudinal- 
transverse splitting.ES 
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FIG. 3: (Color online). Quasiparticle modes of a spinor con- 
densate formed by indirect excitons for negative (a, b) and 
positive (c, d) exchange interaction between excitons. The 
solid lines in the plots (a) and (c) correspond to the disper- 
sions of the quasiparticles without accounting of the spin-orbit 
interaction and numbers point out the bare modes described 
by the formulae ( 20 1- 1 24 1. The accounting of SOI removes the 



spin degeneracy and increases the number of modes. Concen- 
tration of the particles is taken as n = 10® cm~^ . 
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For the negative exchange interaction accounting of 
Rashba SOI leads to the sphtting of degenerate bogolon 
modes /iwi 3 4, which are hnear in the k — kg region and 
behave as bare SOI modes far from the fco point (Fig. 
[Sjjb), red dashed hnes). The gapped mode fiuj2 (Fig. [3](a), 
red dashed hnes) is only slightly renormalized by spin- 
orbit interaction. One should note that accounting of SOI 
leads to the ground state formed by "ij" two-component 
condensate or four-component condensate, and rules out 
the possibility of the formation of single-component con- 
densate as it was mentioned above. 

In the case of positive exchange interaction both modes 
^1.2,4 and huj2 are renormalized in the k = ko point, 
while for large values of k they approach usual Rashba 
like dispersion for upper and lower mode e± = -E^ ± ak 
(Fig. [Sjjd)). The ground state in this case corresponds 
to four-component condensate only. 

The similar situation occurs if only Dresselhaus inter- 
action term present in the Hamiltonian. However, the 
accounting of both Rashba and Dresselhaus SOI requires 
not trivial ground state definition and we leave it as a 
subject for future research. 



IV. REAL SPACE AND TIME DYNAMICS 

The dynamics of the cold exciton droplet can be stud- 
ied using the set of four Gross-Pitaevskii equations for- 
mulated in the previous section ([7])-([l0|. In the present 
article we do not focus on experiment modelling, but 
rather give a qualitative examples of the possible differ- 
ent types of the dynamics of concentration and spin of 
exciton droplets in real space and time. The initial dis- 
tribution of the excitons in a droplet is modelled by 2 
diameter Gaussian wavepacket with maximal concentra- 
tion of the order order n = 10^ cm~^ . When thermalized 
and being far from the hot center, the droplets reveal the 
physics of a cold boson gas, while an additional indirect 
exciton supply is provided in the central region by exter- 
nal current through the structure.'^ In the present paper 
we account for the finite lifetime of the particles (r = 2 
ns) and study the dynamics of the droplets in stationary 
and non-stationary regimes. 

Let us assume that external optical pump is linear po- 
larized and both -1-1 and —1 bright exciton states are 
created. The kinetic terms in the Gross-Pitaevskii equa- 
tions for ^f+i and 'I'_i cause the diffusion of the parti- 
cles from the center of the spot. The strong dipole-dipole 
indirect exciton interaction leads to the repulsion of par- 
ticles from the high concentration regions and manifests 
itself in formation of concentric propagating rings. While 
the aforementioned terms are present in the case of exci- 
ton polariton condensates and have been widely studied, 
the additional terms leading to the transitions between 
bright and dark states are of great importance for the 
case of the indirect excitons. In the following considera- 
tion we describe the effects coming from all terms of this 
kind separately. 
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(a) (b) 

FIG. 4: (Color online), (a) The bright exciton stationary den- 
sity profile for the linearly polarized two-component (-1-1 and 
— 1) condensate with Dresselhaus SOI present. The spin-orbit 
interaction constants are taken as /3 = 3 ^.eV fj.in (red solid 
line) and /3 = 1.5 neV fim (red dashed line). The grey line 
represents the initial density profile of the spot (magnitude 
is reduced by factor of 3), while blue lines point out the re- 
gions of the effective bright-to-dark conversion, (b) The 2D 
density plot for bright excitons showing the "mexican hat" 
profile — 3 fieV /im). 



Firstly, if only one type of SOI (Rashba or Dresselhaus) 
is present, for cylindrically symmetric pumping spot the 
distribution of the intensity of the photoluminescence in 
the real space governed by the concentration of the bright 
states remains cylindrically symmetric. The processes of 
SOI lead to the conversion of the bright states into dark 
which can be observed in the bright exciton density plots 
(Fig. [4]). One can note from the form of SOI operator 
that the most efficient conversion of bright states 
into dark occurs at the points where the density gradient 
is largest (blue arrows in the Fig. |4]ja)). This yields 
the "mexican hat" profile in the near field distribution of 
photoluminescence which depends on SOI strength (solid 
and dashed lines). 

Besides SOI, another mechanism can lead to bright 
to dark exciton conversion. This is the electron or 
hole exchange interaction term leading to the transitions 
{+1, —1} < — > {+2, —2}. It is described by the last term 
in the Gross-Pitaevskii equations ([7])-([To| with the in- 
teraction constant given hy W = Ve + Vh- In the cal- 
culations we consider a GaAs/AlGaAs/GaAs structure 
with 8nm/Anm/8nm QWs, where the exchange inter- 
action constant equals to Ve — Vh — —1.3 fieV ^iv? . 
One can see that the most efficient dark exciton creation 
process takes places for the highest bright exciton den- 
sity regions. Together with repulsive interaction dynam- 
ics, it leads to the appearance of density modulations 
in the radial direction connected to the rings of bright 
exciton concentration. The Fig. [sjja) shows the differ- 
ence of bright exciton density between the case where 
the exchange term is present, {nl^^gf^^), and case where 

it is zero, {n'^right) ■ '^^^ that the exchange term 
leads to the formation of several dips in the bright ex- 
citon density and causes overall flattening of the density 
profile. For the high density of trapped indirect excitons 
bright-to-dark conversion could play important role caus- 
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FIG. 5: (Color online), (a) The effect of bright-dark exciton 
exchange mixing on real space propagation of cold exciton 
droplet. The plot shows the difference between distributions 
of the photoluminescence in cases where bright-dark exciton 
exchange is present (nj^j^^j) and absent (n^^jg^j). (b) The 
density profile of bright excitons with accounting of exchange 
bright-to-dark conversion (dashed line) and without account- 
ing of exchange interaction (solid line). The exchange matrix 
element is equal to Ve ~ —Vdir/5, which corresponds to the 
large separation between QWs. 



ing the gap in the center of photoluminescence profile.'^ 
However, in the typically studied structures with indi- 
rect excitons where exchange interaction is an order of 
magnitude smaller than the direct interaction the density 
modulation is weak and can be expected to be concealed 
by other factors. 

Next, we study the situation when both Rashba and 
Dresselhaus SOI terms present in the Gross-Pitaevskii 
equations, in this case the dispersions of non- interacting 
particles are anisotropic in k-space. 



£±(k) = 



2m 



± ky^a^ + + 2a;9sin(26lk), (25) 



where 0k denotes the angle between the wave vector k 
and X axis. This can lead to the breaking of the cylin- 
drical symmetry in the system as it is illustrated in the 
Fig. [6] The density profile in these plots is governed 
by the interplay of Rashba and Dresselhaus terms, which 
make it cylindrically non-symmetric. The correspond- 
ing evolution of the circular polarization degree shows a 
four-leaf pattern formation. A similar phenomenon was 
observed for spots of cold exciton condensates .'^Sl One 
can note that in the case of equal strength of Rashba 
and Dresselhaus SOI the integrated circular polarization 
is constantly zero. The spin-orbital interaction strongly 
depends on the material and geometry of the sample. 
For the particular 8nm / Anm / 8nm GaAs/AlGaAs het- 
erostructure the constant of Dresselhaus aSOI can be es- 
timated as /3 = 3 fieV /im, while Rashba constant can 
be tuned in wide diapason by the external gate voltage 

Y [M37 38 

Finally, we consider the case of a circularly polarized 
condensate where only one component is pumped. In 
this case the symmetry between +1 — > +2 and — 1 — > — 2 
conversion is removed and it leads to the appearance of 
effective magnetic field in z-direction acting on a bright 




A 4 

x,|im 



FIG. 6: (Color online). Dynamics of a cold exciton droplet 
with both Rashba and Dresselhaus interactions present (/3 = 3 
fieV fim, a = 0.9/3). The set of density plots shows the 
bright exciton (a) and dark exciton (b) density, and the bright 
circular polarization degree determined as p = d^'+i]^ — 



i*-ir)/(i*- 



-ip + i*-ip) 



(c). 




FIG. 7: (Color online). Dynamics of bright (a) and dark (b) 
components under the circularly polarized initial conditions. 
The evolution of exciton density reveals the rotation in time 
of the velocity field for both species, which can be described as 
resulting from the appearance of the effective magnetic field 
provided by the interplay between SOI and . 



and dark components and arising from spin- dependent 
exciton-exciton interactions. Its interplay with SOI of 
Rashba and Dresslhaus types leads to the rotation of the 
densities of bright and dark components in real space 
and time as it is shown at Fig. [7) However, going to the 
stationary regime the density profile coincides with the 
situation for a linear pump (Fig. [Tjjia)). 



V. CONCLUSIONS 

In conclusion, we analyzed the ground state proper- 
ties and spin dynamics in the system of cold indirect 
excitons, accounting for Rashba and Dresselhaus SOI 
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and spin-dependent exciton-exciton interactions. We 
demonstrated that ground state of the cold exciton 
gas is highly degenerate and calculated the dispersions 
of the Bogoliubov excitations in the system. We have 
shown that accounting of spin structure qualitatively 
changes the dynamics of the exciton droplets in the real 
space and can lead to the formation of the rings and 
"mexican hat" structures in spatial distribution of near 
field photoluminescence provided by bright excitons. For 
the case of both Rashba and Dresselhaus terms present 
we have shown the four-leave pattern formation for the 
distribution of circular polarization in the real space. 

Acknowledgement. Wc thank Yuri G. Rubo and A.V. 
Kavokin for the valuable discussions. This work was sup- 
ported by FP7 IRSES projects "SPINMET" and "PO- 
LAPHEN" and Rannis "Center of Excellence in Polari- 
tonics". O.K. acknowledges the help from Eimskip Fund. 
E.B.M. thanks Universidad Autonoma de Mexico for hos- 
pitality. 



where 2x2 blocks Tq read 



2m, ^ e 



2m. ^ h 



a{d%+idl 



a{~dl + idl) \ 



with rue and ruh being electron and hole mass, respec- 
tively, and we denoted the operators acting on electron 
and a hole by indices e and h. The operator T acts on 

the spinor wave function — , \E' ), 

where indices for each component describe the sign 
of spin projection on z axis for electron and a hole. 
Therefore, non-diagonal terms of matrix Tq couple ^'+4. 

and 5* ^ components and are responsible for the spin 

flip transitions between bright and dark excitonic states 
±1 — > ±2. For the exciton it results into bright-to-dark 
states conversion. The operator Tq can be rewritten in 
exciton center-of-mass coordinate frame using transfor- 
mations 



Appendix A: Spin-orbit interaction for indirect 
exciton 

We consider the indirect exciton - composite boson 
consisting of electron and hole in separated QWs which 
are bounded with attractive Coulomb interaction. The 
generic Hamiltonian of indirect exciton with accounting 
of SOI for electron (for instance of Rashba type) can be 
written as 



R = X^e + (1 - X)r/i, 
r = re - r/j. 



where x — ™e/(™e + frih) is the ratio electron and exci- 
ton. Thus, the expression (A3) reads 



H = H„ + V„ 



(Al) 



where Vi„t — — e /47reeov-^^ + ~ r/^p denotes the in- 
teraction between electron and a hole. The kinetic part 
of the Hamiltonian now includes the SOI term for an 
electron and can be represented as 4 x 4 matrix 



To 
to 



(A2) 



with M = rrie + rrih being exciton mass and fi — 
reduced mass. Let us for simplicity account only for first 

matrix Tq acting on the first pair of states and ^ ^. 

The direct multiplication yields 



1^ V2,vl,+^(R, r) + xa{-dx + zay)*_+(R, r) + [-^ + F(r)]vI/++(R, r) + a{-d, + zO,)vI/_+(R, r) = E^++{B., r), 
^V^vl,_+(R,r) + +*5y)*++(R,r) + [-^V^ + y(r)]*_+(R, r) + a(9, + *9,)vI/++(R, r) = i?vI/_+(R,r). 



Substituting the wave functions in the plane wave form 
«'++(R,r) = e*K^0++(r) and*_+(R,r) = e'^^(j>-+{r) 
one can see that the Rashba SOI affects the center-of- 
mass motion of the indirect exciton. Thereby, the ac- 
counting of spin-orbit interaction for electron leads to 
the appearance of SOI acting on the indirect exciton 
reduced by a factor x = We/M denoting the electron- 



exciton mass ratio. The influence of SOI on the relative 
electron-hole motion can be studied within the pertur- 
bation theory and one can show that for the Is state 
of exciton the first order correction is zero. Therefore, 
in the first order approximation one can neglect relative 
motion terms. 
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Appendix B: Interaction matrix elements for 
indirect excitons 



The problem of calculation of indirect exciton inter- 
action constant is reminiscent to the same problem for 
direct excitons. However, one should account that in the 
former case electron and hole layers are separated by cer- 
tain distance and the wave functions of indirect exciton 
are modified. The indirect exciton interaction matrix el- 



ements were calculated in Ref. [35]. However, in the 
following calculations we will follow the approach used 
in Ref. [30] for direct excitons to write integrals in the 
similar form. As was stated before, the spin-dependent 
exciton interactions are based on four types of Feynman 
diagrams: direct interaction, exciton exchange, electron 
exchange and hole exchange. These interactions can be 
written in general form 



Ki^r(Q,Q',q) = 
Kl,(Q,Q',q) = 
K;c/,(Q,Q',q) = 
Kt,,(Q,Q',q) = 



J d^red^rhd^re'd^rh'*Q(re, rh)*Q/ (re, rh')V/(re, rh, rg', rh')*Q+q(re, rh)*Q'-q(re', rh'), 



(Bl) 



J rf^rerf^rhrf^re'd^rh'*Q(re, rh)*Q.(re, rh')V/(re, rh, rg', rh')^'Q+q(re' , rh')*Q'-q(re, rh), 



(B2) 



J d2rerf^rhd^re'd2rh'*Q(re, rh)*Q,(re, rh')V/(re, rh, rg', rh')^'Q+q(re' , rh)*Q'-q(re, rh'), 



(B3) 



J (fred^rh(fre,d^rh,^*(^{re, rh)*Q,(re, rh')V/(re, rh, re', rh')*Q+q(re, rh')*Q'-q(re' , rh). 



(B4) 



where the Coulomb interaction between the electrons and 
holes of different excitons are 

F/(re, rh, re', rh') = y(|re - re'l) + l/(|rh - rh'|)- 
-l/(|re-rh'|)-y(|re'-rh|), 

where minus sign for last two terms account for attractive 
interaction between an electron and a hole. 

The Coulomb interaction for indirect excitons can be 
written as 



yf(re,rh,re',rh') 



47reen 



1 



re - re 



|rh - rh' 



\/(re-rh')2 + L2 V(i-h-re')2 + L2 



where L is a separation distance between centers of cou- 
pled QWs and we used narrow QW approximation. 

There are several possibilities to construct the wave 
function of indirect exciton. As in the usual case of direct 
excitons we consider the 2D motion of bounded electron- 
hole pair but with fixed separation in z direction equal 
to L. In this case it is convenient to separate the exciton 
center of mass motion and relative motion. The general 
form of wave function is 



/(rihR|i) = 



■iKii R|| 



A 



(ril), 



where new coordinates are r|| = r|— rj^ describing relative 
motion and R|| = /3er| -I- f3hJ^'^^ for exciton center-of-mass 
motion. Here {3^^ me/ (rrie + mh), Ph = mh/{me + ruh) 
and A denotes the area of the sample. While the center 
of mass motion is described by the plane wave function, 
the relative motion part of wave function (/)(r|| ) can be 
represented in several different form^^Sl 



as 



02(r|i) 



1 



^/2^^b{b + ro 



: exp 



2b 



(B5) 



, (B6) 



where qb and 2b are quantities associated to indirect ex- 
citon Bohr radii obtained by the variational procedure 
and tq is variational parameter reminiscent to the sepa- 
ration distance between QWs. 

Direct and exciton exchange interaction of indirect ex- 
citons. To calculate the direct dipole-dipole interaction 
of two indirect excitons we choose the second represen- 
tation of the wave function with 0(r||) = 02 (ry). In the 
following derivation we will omit sign || meaning every- 
where the 2D motion. Thus, the integral can be written 
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Vd^r{Q,Q',(^) = 



exp(2ro/6) 



47reeoA2 (27r6(6 + ro))^ 
exp[-iQ(/3ere + Phrn)] exp[-iQ'(^ere' + Ph^h 



d2red^r/id2r^,c?2^^, exp - 

111 1 



) exp 



X exp[«(Q + q)(/3ere + ^S/ir;,)] exp[i(Q' - q)(/3ere' + 



(B7) 



r 



This integral can be simplified if one introduces center- 
of-motion coordinates for both indirect excitons: R — 



PeVe + PhYh, R' = PeYe' + Ph^h' , P = - Th and p' 

Te' - Vh' ■ Then Eq. (IBTI) yields 



Vdir(q) 



exp(2ro/6) 



47reeoA2 (27r6(6 + ro))2 
1 



(fp(fp'(fll(fli' exp 
1 



I 1 J^^Pl 



exp[iq(R- R')]x 



1 



|/3;,(p-p') + R-R'l |-/3e(p-p')+R-R'l V(/3/«P + /9eP' + R-R')' + i' 
1 



v/(-/3eP-/9hP' + R-R')' + i' 



(B8) 



r 



where one can note that complex exponents with Q and 
Q' cancel each other. It is convenient to use the next 



substitutions ^ = R - R', cr (R + R')/2. The integral 
( B8 ) rewritten in a new variables reads 



^r I \ exp(2ro/6) [ ^ ,2 i,2t- ( ^ + ^1 

1 1 1 



) exp 



p " + rl 
b 



exp[iq^] X 
1 



\Mp-p')+i\ |-/3e(p-p')+^l ^{PhP + PeP'+C^+L^ ^{-~P,p-PhP' +$.Y+L-^ 



(B9) 



where the two dimensional integral over variable cr gives 
the area A. The expression ( B9 1 represents a sum of four 



integrals for electron-electron, hole-hole and electron-hole 
mixed interaction 



(BIO) 



where we defined the constant C — -7-^ — j /o^F/^I''^nn2 ■ 



Let us calculate all integrals separately. The first inte- 
gral lee' yields 
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exp 



C / d^pd^p'expl - 



) cxp 



^f5iip-p')^+e 



2ttC / d^p(fp'cxp - 



) cxp 



\/a^ + 7o 



)exp(ig(^ cos( 







(Bll) 



I dVrfV exp ( - V^^itll) exp ( - ^^'^^'° ) exp[-^,qp] exp[/3,qp'], 



where we made substitution = I3^{p — p')^. The inte- 
grals on p and p' are identical and can be factorized into 
Jee'(q) =^p{(l), where 



^p{<i) ^ J dppexp 



71. /-Zti- 
0^ / g-/3hgp 



dpp exp 



Jo(/3,,gp). (B12) 



It is not possible to calculate integral (B12) analytically 



in general case, but we are interested in q = limit. 
Then, the Bessel function of zero order Jq — \ and we 
can change the variable to — ^ r\ and integrate by 
parts 



Ip(q — > 0) = 27r / dppexp ( 



b 



27r 



+ 00 



Finally, the integral lee' reads 

I„c.-o ^ ^ _exp(2ro/b)^ 2^^^^^,,„/,^^^_^^^^^. 



27reeoA(7 



One can see that expression for integral Ihh' coincides 
with lee' with substitution f3fi — ^ /3e- Therefore, in the 
q — > limit they are equal Ihh' — /2TreeoAq. 



Now let us calculate the second type of integrals re- 
sponsible for the electron-hole attractive interaction 



Ie/.'(q) = C J d^pd^p'd^^exp ( - 

i 



X exp I 



p'' + rl 



cxp[iq|] 



v/(/?,.p + /3eP'+^)2+L2- 



and performing the substitution x = PhP + PeP' + ^ we 
can write 



^eh' (q) = C y d^pd^p'd^x exp ( ^ 
C / d'pd^p' exp ( - ^^lill ) exp ( 



exp 



exp[-iqx] exp[-i/3/iqp] exp[-i/3eqp'] 



1 



cxp[-i/3;,qp]exp[-i/3eqp'] / dx . = / 

JQ VX +i "'0 



2ir 



c 



27re-9^ 
9 



/ £-pSf> exp !i j exp y - -^^ j exp[-i/3,,qp] exp[-i/3eqp']. 



/,g«9X cos 



(B13) 
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The same factorization can be done as in the integral 



(Bll) and expression (B13) yields 



leh' (q ^ 0) = 



e2 e-?^ 



27reeo^ 



(B14) 



The similar result was obtained using another approach 
by authors in Ref. [32] , 



One can check that the same expression is valid for 
^/ie'(q — > 0). Finally, the sum of four integrals yields 



L. 



(B15) 



The exciton exchange interaction V^^^ can be written 
in the same way as direct interaction from the general 
form (|B2| 



Kl,(Q,Q',q) 



exp(2ro/6) 



exp[-iQ(/3ere + hrh)] exp[-iQ'(/3ere' + PhYw)] 



d''Y^d^Vh(I''Ye'd''Yh, exp 
1 1 



exp 



1 



exp[i(Q + q)(/3ere' + h^h')] exp[i(Q' - q)(/3ere + PhVh 

( 



exp 



cxp[-iAQ(R - R')] exp[-iq(R - R')] 

1 



exp(2ro/b) 
47reeoA2 {2Tih{h + ro)Y 

1 



d'^pd'^p'd^Rd^R' exp ( 
1 



V(/3/>P + /3eP' + R-R')'+i' V(-/?eP - PhP' + R - R')' + 



where we defined exchanged momentum between elec- 
trons as AQ — Q — Q'. One can see that for small 
exchanged momentum between exciton which is true for 
weakly interacting exciton gas the expression for V^^f^ 
coincides with Vdir for g — > 0. 



Electron and hole exchange interaction of indirect ex- 
citons. The general form of electron exchange interaction 
for indirect excitons is given by Eq. ( B3 ) . For the calcu- 



lations it is more convenient to choose the indirect exci- 
ton wave function in the form ( |B5[ ). Thus the exchange 
interaction matrix element yields 



Kxch(Q,Q',q) 



d^Ted^rhd^re'd^rhr exp 



a-B 



exp 



a-B 



X exp 



a-B 



-) exp 



gb 



exp[j(Q + q)(/3ere + Phrh)] cxp[z(Q' - qjiPeVe' + l^hYh' 



exp[-iQ(/3ere -f Phr^h)] exp[-iQ'(^ei"e' + PhYh')]^- 
1 1 1 



The exact calculation of exchange integral is straightfor- 
ward but tedious. Using the same steps as for direct in- 
teraction calculation and performing in the end substitu- 
tions yi ^ {^-l3ep~l3hp')laB, yi = {^+PhP+l3eP')/aB, 
X = p/gb and L = L/aB one gets the final expression of 



electron exchange interaction 

K'Lch = -^^^ «B-Ie'.e/.(AQ,<Z,e,/?e), (B16) 

where the exchange integral is given by 
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f27r />C30 j>2'K />C30 j>2'K 

dx I dQx / dyi dQi dy2 / ^62 ■ 2:2/12/2 cos 























AQaB[l3eX cos{Q - 0^) + 



PeVi cos(e - 61)] + gasf-xcose^ - P^yi cos 61 + (1 - /3e2/2 cos 62)] exp(-Y + L2)exp(- J yl + 



exp(— a/?/| + L"^) exp(— Y (2/2 cos 62 — 2/1 cos 9i — xcos Qx)^ + (2/2 sin 62 — 2/1 sin 61 — xsinQxY + L'^)"- 



^/yf + x^ + 22/ia; cos(ei - 6^;) ^2/2 + + 22/22; cos(e2 - 6:^) ^^2 ^ ^2 _^ 



L2 



(B17) 



r 



with 9 being the angle between AQ and q. We are 
interested in the case when AQ = and g = 0. The cal- 
culation of exchange integral thus requires numerical in- 
tegration with multidimensional Monte Carlo algorithm. 



Moreover, it is obvious that similarly to the case of direct 
excitons electron and holes exchange interactions have 
the same value for g — > 0. 



^ L. Pitaevskii and S. Stringari, Bose-Einstein Condensation 
(Oxford University Press, Oxford, 2003) and A. Legett, 
Quantum Liquids: Base Condensation and Cooper Pairing 
in Condensed-Matter Systems (Oxford University Press, 
Oxford, 2006). 

^ S. A. Moskalenko and D. W. Snoke, Bose-Einstein Conden- 
sation of Excitons and Biexcitons (Cambridge University 
Press, Cambridge, U.K., 2000). 

^ L. V. Keldysh and Yu. V. Kopaev, Sov. Phys. Solid State 
6, 2219 (1965); L. V. Keldysh and A. N. Kozlov, Sov. Phys. 
JETP 27, 521 (1968). 

* I. B. Spielman, J. P. Eisenstein, L. N. Pfeiffer, and K. W. 
West, Phys. Rev. Lett. 87, 036803 (2001); J.P. Eisenstein 
and A. H. MacDonald, Nature 432, 691 (2004). 

^ S. O. Demokritov, V. E. Demidov, O. Dzyapko, G. A. 
Melkov, A. A. Serga, B. Hillebrands, A. N. Slavin, Nature 
443, 430 (2006). 

® J. Kasprzak, M. Richard, S. Kundermann, A. Baas, P. 
Jeambrun, J. M. J. Keeling, F. M. Marchetti, M. H. Szy- 
manska, R. Andre, J. L. Staehli, V. Savona, P. B. Lit- 
tlewood, B. Deveaud, and Le Si Dang, Nature 443 , 409 
(2006). 

^ R. Balili, V. Hartwell, D. Snoke, L. Pfeiffer, and K. West, 

Science 316, 1007 (2007). 
® J. J. Baumberg, A. V. Kavokin, S. Christopoulos, A. J. 

D. Grundy, R. Butte, G. Christmann, D. D. Solnyshkov, 

G. Malpuech, G. Baldassarri Hoger von Hogersthal, E. 

Feltin, J.-F. Carlin, and N. Grandjean, Phys. Rev. Lett. 

101, 136409 (2008). 
® V. B. Timofeev, Physics-Uspekhi 48, 295 (2005). 
1° L. V. Butov, J. Phys.: Condens. Matter 19, 295202 (2007). 
" L. V. Butov et al., Phys. Rev. Lett. 86, 5608 (2001). 
^2 D. Snoke, Science 298, 1368 (2002). 

Yu. E. Lozovik and V. I. Yudson, Zh. Eksp. Teor. Fiz. 

Pisma Red. 22, 26 (1975) [JETP Lett. 22, 26 (1975)]; Zh. 

Eksp. Teor. Fiz. 71, 738 (1976) [Sov. Phys. JETP 44, 389 

(1976)]; Solid State Commun. 18, 628 (1976). 
" S. I. Shevchenko, Phys. Rev. Lett. 72, 3242 (1994). 

X. Zhu, P. B. Littlewood, M. S. Hybertsen, and T. M. Rice, 

Phys. Rev. Lett. 74, 1633 (1995). 



O. L. Berman, Y. E. Lozovik, D. W. Snoke, and R. D. 
Coalson, Phys. Rev. B 70, 235310 (2004). 
^'^ V. Larionov, V. B. Timofeev, J. Hvam, and K. Soerensen, 
Zh. Eksp. Teor. Fiz. 117, 1255 (2000) [JETP 90, 1093 

(2000) ]; A. V. Larionov and V. B. Timofeev, Pisma Zh. 
Eksp. Teor. Fiz. 73, 342 (2001) [JETP Lett. 73, 301 

(2001) ]. 

1* D. Snoke, S. Denev, Y. Liu, L. Pfeiffer, and K. West, Na- 
ture 418, 754 (2002). 

^® L. V. Butov, C. W. Lai, A. L. Ivanov, A. C. Gossard, and 
D. S. Chemla, Nature 417, 47 (2002); L. V. Butov, A. C. 
Gossard, and D. S. Chemla, Nature 418, 751 (2002). 
For review of spin dynamics of cavity polaritons see I. A. 
Shelykh, A. V. Kavokin, Yu. G. Rubo, T. C. H. Liew and 
G. Malpuech, Semicond. Sci. Technol. 25, 013001 (2010). 
T. C. H. Liew, I. A. Shelykh, and G. Malpuech, Polaritonic 
Devices, Physica E 43, 1543 (2011). 

22 N. A. Gippius, I. A. Shelykh, D. D. Solnyshkov, S. S. 
Gavrilov, Yu. G. Rubo, A. V. Kavokin, S. G. Tikhodeev, 
G. Malpuech, Phys. Rev. Lett. 98, 236401 (2007). 

2^ Yu. G. Rubo, A. V. Kavokin, I. A. Shelykh, Phys. Lett. A 
358, 227 (2006). 

2" O. Kyriienko, I. A. Shelykh, Phys. Rev. B 84, 125313 
(2011). 

Yu. G. Rubo, A. V. Kavokin, Phys. Rev. B 84, 045309 

(2011). 

A. A. High et al, arXiv:1103.0321 

M. Z. Maialle, E. A. Andrada e Silva, L. J. Sham, Phys. 
Rev. B 47, 15776 (1993). 

R. Winkler, Spin- Orbit Coupling Effects in Two- 

Dimensional Electron and Hole Systems, Springer Tracts 

in Modern Physics, 191, (2003). 
2^ R. H. Silsbee, J. Phys. Condens. Matter 16, R179 (2004). 
^° M. Combescot, O. Betbeder-Matibet, F. Dubin, Physics 

Reports 463, 215 (2008). 

M. Combescot, O. Betbeder-Matibet, and R. Combescot, 
Phys. Rev. B 75, 174305 (2007). 

S. Ben-Tabou de-Leon and B. Laikhtman, Phys. Rev. B 
63, 125306 (2001). 

S. de-Leon and B. Laikhtman, Phys. Rev. B 61, 2874 



14 



(2000). 

E. M. Lifshitz and L. P. Pitaevskii, Statistical Physics, part 

2 (Pergamon Press, New York, 1980). 

I. A. Shelykh, Yuri G. Rubo, G. Malpuech, D. D. Sol- 

nyshkov, and A. Kavokin, Phys. Rev. Lett. 97, 066402 

(2006). 

A. A. High, J. R. Leonard, A. T. Hammack, M. M. Fogler, 
L. V. Butov, A. V. Kavokin, K. L. Campman, and A. C. 



Gossard, arXiv:1109.0253v2. 

E. Ya. Sherman, A. Najmaie, and J. E. Sipe, Appl. Phys. 
Lett. 86, 122103 (2005). 

S. Giglberger et al., Phys. Rev. B 75, 035327 (2007). 
N. W. Sinclair et al, Phys. Rev. B 83, 245304 (2011). 
C. Ciuti, V. Savona, C. Piemarocchi, and A. Quattropani, 
Phys. Rev. B 58, 7926 (1998). 



